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ABSTRACT 



We have studied the fascinating dynamics of the nearby Vela pulsar's nebula in a 
campaign comprising eleven 40 ks observations with Chandra X-ray Observatory ( CXO) . 
The deepest yet images revealed the shape, structure, and motion of the 2-arcminute- 
long pulsar jet. We find that the jet's shape and dynamics are remarkably consistent 
with that of a steadily turning helix projected on the sky. We discuss possible implica- 
tions of our results, including free precession of the neutron star and MHD instability 
scenarios. 

Subject headings: pulsars: individual (Vela); outflows 



1. Introduction 



The rapid rotation of a pulsar's magnetosphere gives rise to strong electric fields, particle 
acceleration and pulsed radiation across the whole EM spectrum. Pulsar wind and associated radi- 
ation can fill a substantial volume in the vicinity of a pulsar and be observ able throughout the EM 
spectrum as a pulsar wind nebula (PWN; for recent reviews of PWNe see lGaensler &: Stand 120061 ; 
Kargaltsev Pavlov! 120081 ) . The termination shock, where bulk flow kinetic energy is converted 
into internal energy of the plasma, is resolved for th e case of the Crab Nebula (and some other 



PWN e) as the inner edge of the torus seen in X-rays (jWeisskopf et al 



2000 



Kargaltsev &: Pavlov 



2008 ). In addition to toroidal stru ctures, many PWNe exhibi t jets which are believed to be out- 
flows along the pulsar's spin axis (jKargaltsev Pavlovl 120081 ) . The formation, confinement and 
structure of jets are still poorly understood, particularly for pulsar jet s, which are very different 
from AGN, X-ray binary and protostellar jets in some respects (see, e.g., lHardedl2008l for a review). 
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The proximity of the Vela PWN offers a unique opportunity to carry out a detailed, time-resolved 
study of jet dynamics on timescales of months. Such studies can shed new light on jet physics in 
general, but they are impossible to perform for AGNs jets, which are the most commonly discussed 
examples of magnetized MHD outflows. 

The Vela pulsar is a young (pulsation period P = 89ms and spin-down age r = 11.3 kyr), 
energetic pulsar locat ed at a distance d ks 300 p c. The pulsar is fam ous for its 7-ray brightness 
(jHartman et al.lll999l ) and glitching behavior (IMcCulloch et al.lll990l). The lat ter was considered 
strong evidence of superfluidity in the pulsar core (jSedrakian fc Corded 1 19991). The Vela pulsar 



powers the nearest bright X-ray PWN, with a puzzling double-arc structure (jHelfand et all 12001 
Pavlov et al]l2001bll20oi ). 



Previous X-ray observations with the Chandra X-ray Observatory ( CXO ) led to the discovery 
of the pulsar's outer jet (beyond the arcs, in the direction of the pro per motion), which i s fain t 
compared to the inner nebula, and has a peculiar curved morphology (jPavlov et al.l l2001al . 120031 ) . 
Unfortunately, the read-out streak, the chip gaps, and the suboptimal cadence hampered modeling 
of the outflow dynamics. Here we present more recent CXO follow-up observations (3 images in 
2009 and 8 in 2010, a total of ~440 ks of new data), with each image twice as deep as in the previous 
observations. In this Letter we report the first results of these observations. 



2. Observations and Data reduction 

We observed the Vela Nebula with the Advanced CCD Imaging Spectrometer (ACIS) aboard 
Chandra. Three observations were performed in July 2009 and eight in June-September 2010. Each 
observation had an exposure time of about 40 ks, and the intervals between successive observations 
were about a week. The FAINT telemetry mode, with a frame time of 3.24 s, was used with the 
target imaged on the back-illuminated S3 chip. The read-out streak (caused by the bright pulsar) 
did not overlap with the jet. We applied standard filtering to the recorded events and ignored 
energies outside the 0.5-8 keV band. The sequence of the resultant images is shown in Figure [T]l2j 



3. Analysis and Model 



The jet in the direction of the pulsar proper motion (northwest) has two components: the 
bright inner jet, which extends for ~ 13" from the pulsar up to the outer arc of the PWN, and the 
outer jet, which is much fainter but extends further for another ~ 100". For brevity, hereafter we 
will refer to the outer northwe st (NW) jet as "the jet" unless it causes ambiguity. With new deeper 
images (compared to those in lPavlov et al.l 120031 ) . the fainter emission filling the jet is clearly seen 
in addition to the brighter knots (see Figures [T] and [2]) . Despite the changing shape, the jet as a 
whole remains confined to a fairly well defined region (jet "envelope"), which is initially narrow and 
pointing in the direction of the pulsar's proper motion (as does the inner jet), but then gradually 
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Fig. 1. — The 2009 data with the projected trace of the best-fit helical model superimposed (green). 
The envelope and guiding line (see text) are also marked (dotted and dashed blue lines, respec- 
tively). A black X marks the location of the pulsar and a white vector the direction of proper 
motion on each image. The parameters for the helix are r = 122 day, A = 0.23 ly and an opening 
half-angle of 5°. 




Fig. 2.— As FigureH but for the 2010 data. 
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bends towards the southwest (SW) at greater distances from the pulsar (see Figures 1 and 2, blue 
curves; also below). 

The jet shape and dynamics (see Figures 1 and 2 and the corresponding animation) strikingly 
resemble that of a rotating helix projected on the sky. The successive positions of the knots detected 
within the jet in several images, do not appear to have a large velocity component transverse to 
the jet but rather move approximately radially away from the pulsar. The brighter inner jet does 
not show nearly as much motion as the outer jet. H owever, the point at which the outer jet 



emerges beyond the outer arc (see lPavlov et al.ll2001bl for definitions) changes position with time 



moving back and forth along the outer arc edge. This remarkable phenomenology prompted us to 
explore a simple geometrical model where a small perturbation, which initially occurs close to the 
base of the jet, propagates ballistically and grows with the distance from the pulsar resulting in a 
three-dimensional helix, which we view projected onto the sky. 



In our model we adopt the geometric orientation suggested by iHelfand et al.l (|200ll ). which 
implies that the angle between the axis of the inner jet (the latter is directed along the pulsar's 
spin axis and the proper motion direction) and the line-of-sight (LoS, normal to the sky) is 53°, 
with the inner jet pointed away from the observer (i.e., behind the plane of the sky). The ballistic 
motion of the material in the (o uter) jet is gradually modified due to the global bending of the jet 



(discussed by lPavlov et al.112003 ; see also the dashed blue curve in Figures 1 and 2). We require the 
helix to be confined within a bent envelope (see blue curves in Figures 1 and 2), which we define 
at the outer boundary of the jet emission on the sky, from the sum of all 11 images. We define 
the "guiding" line as a curve going through the middle of the envelope. The helix radius increases 
with distance from the pulsar (forming a cone with opening half-angle 5° ± 1°, as measured from 
the pulsar). The helix and its trace (projection) on the sky were constructed in the following way. 
The helix was defined parametrically using a rotating vector and a moving circle. At each point 
along the guiding line (dashed blue curve in Figures 1 and 2) we defined a circle in a plane locally 
orthogonal to the guiding line. The circle radii are increasing with distance from the pulsar, so that 
each of the circles meets the jet envelope (dotted curves in Figures 1 and 2). For each circle we 
drew a vector from the circle's center to a point on the circle. We then rotated these vectors in their 
circle planes so that the azimuthal angle (in the circle plane) changes linearly with the distance 
from the pulsar. The endpoints of these vectors trace a 3-dimensional helix, which we projected 
onto the sky. The dynamics was introduced by having the azimuthal angle of the first vector (near 
the pulsar) change linearly with time, completing one turn in time r. Thus, by construction, in our 
model the speed of a helical pattern is constant, v = A/r, where A is the helix pitch (the distance 
between successive turns of the helix). 

As a next step, we fitted our dynamical model to the CXO image series by varying r, A, and 
the initial azimuthal angle of the vector near the pulsar (which sets the jet "launching" direction 
at a given time). The model parameters were optimized to maximize the number of photons 
falling in a 3pix-wide stripe along the trace. The best match to the data corresponds to the 
stripe containing ~ 500 extra counts (~4200 counts total) compared to the average value for 
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randomly chosen parameters (~ 3700 counts total). To obtain the most conservative estimate of 
this significance we performed extensive Monte-Carlo simulations. In particular we simulated the 
jet shapes for 50< r <200d and 0.1< A <0.31y and any initial azimuthal angle, which are also 
the ranges we used in our best-fit search. The simulations resulted in 10000 realizations for the jet 
shape. We then calculated the ratio of the number of realizations where the jet's trace intercepts 
>4200 photons to the total number of MC realizations. This resulted in the chance coincidence 
probability of 0.4% which corresponds to ~ 3a significance. A simplistic estimate based on Poisson 
counting statistics leads to larger significance. Because the data consist of two series with about 
a one-year gap in-between, there is strong aliasing in the goodness-of-fit measure as a function of 
t, resulting in several acceptable periods, r = 73 ± 2,91 ± 5, 129 ± 10 d. We repeated this analysis 
including older observations from 2000-2002. The Monte-Carlo significance did not improve, but 
one of the above solutions was preferred above the others: A = 0.230 =L 0.014 ly (corresponding 
to a projected angular distance of 40" on the sky) and r = 122 ± 5d (~3 cycles/year), implying 
v = A/t = 0.23cyr/0.33 yr = 0.7c. This flow speed is in agreemen t with the earlier e stimates of 
the apparent speed of individual knots in the jet, v = 0.3c — 0.7c (IPavlov et al.ll2003l ). We have 
sketched the helical trace corresponding to the best fit on the images in Figures [T] and [2] (green 
curves) . 



4. Discussion 



The model described in Section 3 appears to provide a simple description of the observed jet 
geometry and dynamics. According to our model, the flow of particles is ballistic: the bulk velocity 
is predominantly radial (away from the pulsar) in the direction defined by its launching angle, 
and the helix arises only from the steady change of that launching direction (similar to a rotary 
garden sprinkler). Indeed, the individual knots within the jet appear to move away from the pulsar 
radially, rather than following the helical shape of the jet (Figures 1 and 2). Below we discuss the 
implications of a helical structure produced by a ballistic flow. 

Helical patterns are often seen in AGN j ets (e.g.JPerucho et al ] l2012h . protostars (e.g-. lLee et al. 



20071 ) and accreting binary systems (SS443; IStirling et al.ll2002l ). Among the possible mechanisms 
leading to a helical structure, the most frequentl y mentioned are the helical instability mode in a 
magnetized plasma column (the 'kink' i nstability; Harde ell201ll . 120061 and references therein) or jet 



2011J 

precession (|Vernaleo fc Reynolds] 12006] ; iMasciadri et al.ll2002l ). Perhaps the clearest example of a 
Galactic helical jet is that of the micro-quasar SS 433, where the shape of the jet has been explained 
as the motion of a matter stream ej ected from a precessing nozzle located in the vicinity of the 
compact object (IStirling et al.ll2002l ). Precession in this latter case is torque-driven through the 
interaction of the compact object with its binary companion. Al so in AGN jet models involving 



precession ([Abraham &; Romero 



1999 



Tatevama &: Kinghamll2004l ) , the precession is torque-driven 



via interactions between the accretor and the acc retion disk. An alternative magneto- hydrodynamic 



(MHD) instability scenario (jHardee 



2011 



20061 ) does not require any external force: rather, the 
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helical pattern grows from a small perturbation which may be any random fluctuation. 

For the jet of the Vela pulsar (which lacks any evidence of accretion or a companion), we 
consider two mechanisms which could give rise to helical structure: free precession (i.e., without 
external torque), and a kink MHD instability in the jet flow. Whilst the first is conceptually simple, 
requiring only oblateness of the pulsar, its existence has not yet been established in any pulsar 
(although several claims have been made, see below). The latter option has undergone extensive 
theoretical and numerical investigations for AGN jets, but not pulsar jets, whose properties are 
different. Moreover, the two mechanisms need not be mutually exclusive because the initial small 
perturbation in the jet flow can be caused by precession and then be amplified via kink instability. 

The mechanism responsible for pulsar jet formation is a matter of debate. If the jets are formed 
at a relatively large distances from a pulsar, e.g., via a magnetically channeled backflow of particles 
from the equatorial post-shock region^], then one might expect the effects of a free precession on 
a jet be somewhat washed out. On the other hand, it is possible that despite t he presence of the 



backflow at larger distances, the jet is actually formed much closer to the pulsar (|Lvubarskvl l2012) 
and then the impact of precession on the jet's behavior would be more pronounced. 



4.1. Helical structure due to precession 



The Vela pulsar spin axis may be freely precessing with the period we measured for the jet 
from our model, in which case the jet is simply launched ballistically along the spin axis. Pre- 
vious hints of free precession in pulsars have been reported, based on (quas i -)per iodic changes in 

). The periodic 



pulsar spectra, pulse profiles or timing residuals (jrlaberl et al 



200e 



signature seen in PSR B 1828— 11 was interpreted as precession (jStairs et al- 



ii; Jonesl 



2011 



2000 ), but after ten 



years of further observation, this interpretation has not been confirmed. The varying X-ray flux 
of the isolated neutron star RX J0720.4— 3125 was suggested to be due to precession but, a gain, 



further observations have failed to firmly establish this conclusion (Ho hle et al. 



2009 



2012); the 



perio d ic flux variati ons could plausibly be caused by effects other than precession (jZhang fc Xie 
2012J). iPopovl (|2008l ). commenting on the putative precession of PSR B1828— 11, argued for a rela- 
tionship between the precession perio d r and the p ulsar spin period P as r ~ kP 1 ^ 2 , which would 
predict r = 120 d for Vela. Similarly, Ijonesl (120111 ) measured the periodic modulations of several 
pulsars, finding a correlation between the long-term period and the pulse period, which could be 
ascribed to precession. Periods in the 80-200 d range are predicted for a pulsar with P=89ms (see 
Figure 1 in |j( 



Jones 2011 



Therefore, it is conceivable that the helical shape of the Vela pulsar jet 
could be due to free prec ession. If the Vela pulsar is ind eed precessing, long-term timing in radio 
may be able to detect it (jDeshpande Sz McCullochlll9 96). but such detections are complicated by 



1 A possibility suggested bv lBogovalov fc Khangoulianl J2002h : iKhangoulian fc Bogovalovl J2OO3I ) and reproduced 
in the numerical simulations of iKomissarov fc Lvubarskvl 1)20041 ): iDel Zanna et alJ (|2004l . I2006T ) . 
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glitching. If precession is confirmed f or the Vela pulsa r, it should be present in other pulsars as 
well. Interestingly, recent analysis by IWeisskopa (|201ll ) shows that the Crab pulsar jet also has a 
corkscrew structure which evolves with time. 

The precession of a pulsar has interesting theoretical implications. For our helical model, we 
can infer the stellar oblateness parameter e = (I3 — = P^m /Pprec ~ 8 x 10~ 9 (where I\ and 

-Z3 are principal moments of inertia). As pointed out by I Jones Anderssonl (|200ll ). this imposes 



a limit on how much of the pulsar mass can be in the superfluid interior. One might expect the 
superfuidity to damp any precession effects, but the cou pling between crust and core is co mplex, 
and th e core may partially participate in the precession (jLevin &: D'Angeloll2004l ). Indeed. I Popov 
(|2008l ) suggests that glitching - exchange of angular momentum between the crust and the core - 
can actually be the excitation mechanism for precession. If so, then it is not surprising that Vela, 
which is such a prolific glitcher, should show the clearest signs of precession. 

We would also like to point out that the rest of the bright inner PWN also appears to be 
variable on a similar timescale. For example, Figure [3] shows the difference image throughout the 
2010 data (last image of the series minus the first image), revealing the changes in the inner arc-like 
structures and the inner south-eastern jet. Changes, labelled in Figure O can be easily seen. We 
defer detailed investigation of the geometric and spectral evolution of the inner nebula for a separate 
paper. It should be noted, however, that the only mechanism capable of prod ucing a coordinated 
(coherent) change in the arcs (likely brightened limbs of the equatorial outflow; iPavlov et al.l l2003) 
and both jets must be acting at the origin of these outflows (i.e. at the pulsar) making the 
precession a natural candidate. Changes in the Vela PWN morphology seen in earlier CXO images 



were suggested to be caused by precession (Deshpande & Radhakrishnan 



earlier (quasi-)periodic timing signatures (jDeshpande fc M cCulloch 



20071 ) , perhaps related to 



1996 ), suggesting a p recession 



period of 330 d. We note that the Vela PWN model of lDeshpande fc Radhakrishnan! (|2007l ) assumes 
that the arcs are "traces of the particle beams from the two magnetic poles at the shock front" and 
the jet is a "manifestation of a physical flow of particles along the magnetic axis of the star, not 
along the rotation axis." Although their assumptions clearly differ from those in our model (where 
the jet is along the pulsar's spin axis), the observational implications of the free precession scenario 
could be similar (regardless of the specific model for the PWN features), such as correlated periodic 
variability in the PWN structure. 



4.2. Helical structure due to kink instability 



Another possible scenario for the observed jet shape and dynamics is the kink instability. 
Although our helical model has strict periodicity built-in, the model may still approximate, to some 
extent, the behavior expected due to the kink depending on the nature of the initial perturbation 
(which can also be the precession) . The M HD helical kink instability is often invoked to interpret the 



AGN jet shapes (e.g.. lHardee 



2011 



20061 and references therein) . Numerical simulations of the kink 



instability in the non-linear regime can produce structures resembling a helix, with current driven, 
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velocity s hear driven, o r jet precession mechanisms being considered as viable options for driving 
the kink (|Hardedl201ll ). Pulsar jets were discovered relatively recently, their physical properties 



are still poorly known, and their modeling is still in a nascent state compared to the AGN jets. 
Although the basic MHD principles should probably be applicable to pulsar jets, existing AGN jet 
models may not be directly applicable, because the properties of these systems (e.g., baryon loading, 
bulk Lorentz factor, magnetic-to-kinetic energy ratio, jet formation mechanism, jet confinement 
mechanism) are likely quite different. Even in the linear approximation the analysis for the local 
i nstability mode s in pulsar jets critically depends on the poorly known magnetization parameter 
(|Begelmanl 11998). Still, some results from the numerical modeling of AGN jets might be useful 
for qualitative understanding of pulsar jets. For the MHD instability, we can crudely associate the 
period in our helical model, r, with the characteristic timescale of the appearance of successive 
kinks (which will not be truly periodi c unless trigge r ed by prec ession) . According to the recent 
MHD simulations done for AGN jets bv lMizuno et al.l (|2009l . l2012l ). the instability growth timescale 
for the global kink development (i.e., when the jet resembles a spiral with a radius comparable to 
the jet radius) can be crudely estimated as U ~ WRj et /vA, where Rj e t is the jet radius and va 
is the Alfven speed. Using t he measured jet ra dius Rj e t ~ 5" and the Alfven speed estimated 
for equipartition, va = 0.77c (jPavlov et al.l 12003) . we find the instability growth timescale to be 
~ 100 ds. In the images shown in Figures 1 and 2 the distinct helical pattern typically emerges 
by a distance dj ~ 1' from the pulsar, which implies the flow velocity di/ti ~ 0.9c assuming that 
the instability starts at much smaller distance from the pulsar than d{. Although the implied 
flow speed is somewhat higher than the flow speed inferred above and the apparent speeds of the 
knots (Pavlov et al. 2003), these are consistent given the uncertainties involved, which support the 
above-described scenario, if the jet's magnetic field strength is near equipartiton. 

We note, however, that the jet retained approximately the same helical morphology between 
the 2010 and 2011 series which corresponds to over 400 d. On this substantially larger timescale, 
one could expect to see highly nonlinear instability development, leading to the jet disruption. 
However, with the data available we cannot tell whether the jet shape remained stable from 2009 
to 2010 or if it had disrupted and regrown. We also note that t he extreme bending was observed 
only once, in the 2000 February image (see iPavlov et al.l 120031 ) . and it might have been due to 
a disruptive, highly-developed kink. The rareness of such events hints at a substantially larger 
instability timescale than the estimate above. In this case the origin of the stable helical pattern 
may be better explained by pure precession. On the other hand, there may be factors leading to the 
delay in the nonlinear instability growth, such as the pressure of the ambient medium becoming 
important far enough from the pulsar. At large distances from the pulsar the jet behavior is 
also affected by the ram pressure due to the combined action of the pulsar's motion and a strong 
wind within the SNR, poss ibly caused by the passage of the reverse shock (Pavlov et al., 2003; 
Chevalier &: Revnoldsll201ll ). This ram pressure can explain the clockwise bending of the jet's end 



and the distortion of the helical structure at large distances from the pulsar. 




Fig. 3. — The difference image between the first and last observation of the 2010 series. We label 
the axis of the inner jets (blue dotted line), the pulsar position (blue circle), and the arc positions 
in the first image (darker region with green dashed line) and last image (brighter region with cyan 
dashed line). Green arrows show the direction of apparent motion, in both the arcs and the jet. 
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5. Conclusion 

The Vela pulsar system continues to be a fascinating testbed for magnetized outflow models. 
It may also become the first NS known to precess freely, if our results are confirmed by further 
modeling of the inner PWN, by additional X-ray observations, and/or by careful pulsar timing 
analysis, which we strongly encourage. If the precession is confirmed, the Vela pulsar may become 
one of the best persistent sources of gravitational waves for the next generation of detectors, given 
that current detecting technology is most sensitive to signal freque ncies >10Hz, and sensitivity is 



greatly enhanced by searching at a known frequency and position (jPitkin et al.1 120111 ). Numerical 
MHD simulations of pulsar jets are also required to make further progress and to reveal the physical 
conditions in the jet and instability development. 
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